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ABSTRACT 

^ ■ Old neutron stars (ONSs) which have radiated away their internal and rotational 

energy may still shine if accreting the interstellar medium. Despite their large number, 
only two promising candidates have been detected so far and rather stringent limits 
on their observability follow from the analysis of ROSAT surveys. This contrasts with 
optimistic theoretical estimates that predicted a large number of sources in ROSAT 
fields. We have reconsidered the issue of ONSs observability, accounting for the spin 
and magnetic field evolution over the neutron star lifetime. In the framework of a 
spin-induced field decay model, we show that the total number of ONSs which are, at 
present, in the accretion stage is reduced by a factor ~ 5 over previous figures if the 
characteristic timescale for crustal current dissipation is ~ 10 8 — 10 9 yr. This brings 
theoretical predictions much closer to observational limits. Most ONSs should be at 
present in the propeller phase and, if subject to episodic flaring, they could be observable. 
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1. Introduction 

Old neutron stars (ONSs), i.e. neutron stars which 
have evolved beyond the pulsar phase, are expected 
to be quite a large Galactic population counting as 
many as N ~ 10 9 objects. Ostriker, Rees, & Silk 
(1970) were the first to suggest that accretion of the 
interstellar medium (ISM) may produce enough lu- 
minosity to make the closest stars observable and 
that ONSs may contribute to the soft X-ray back- 
ground. More recently, Treves, & Colpi (1991) and 
Blacs, & Madau (1993) have shown that several thou- 
sands ONSs should be present in the ROSAT All 
Sky Survey. Rather optimistic predictions for the 
number of detectable sources were also presented 
by Madau, & Blaes (1994) and Colpi, Campana, & 
Treves (1993) for ONSs embedded in Giant Molec- 
ular Clouds (GMC) and by Zane et al. (1996) for 
ONSs in the solar proximity. However, despite the 
intense observational efforts, the search for ONSs pro- 
duced, up to now, just a handful of candidates, out 
of which only two, RX J18653. 5-3754 (Walter, Wolk, 
& Neuhauser 1996) and RX J0720.4-3125 (Haberl et 
al. 1996, 1997), seem indeed promising. Moreover, 
recent analyses of ROSAT fields in the direction of 
GMCs (Belloni, Zampieri, & Campana 1996; Motch et 
al. 1997) placed rather stringent upper limits on the 
number of ONSs in ROSAT images. In all these cases, 
the predicted number of ONSs turns out to be too 
close to the observed number of non optically iden- 
tified sources (NOIDs) which should comprise many 
other types of objects, like white dwarfs, active coro- 
nae and AGNs. In the field investigated by Motch 
et al., for example, ~ 5 sources should be detectable 
above 0.02 counts -1 (Zane et al. 1995) while the num- 
ber of NOIDs, at the same limiting flux, is 8. Since 
these are sites where accretion should be the most ef- 
fective, observations seem to suggest that current the- 
oretical predictions are in excess of at least a factor 
~ 5-10. Moreover, it was pointed out by Zane et al. 
(1995) that most optimistic models produce a source 
density of - 10 deg~ 2 in ROSAT PSPC pointings 
at the limiting flux of 10~ 3 count s _1 , embarrassingly 
close to the average density of NOIDs, ~ 30 deg -2 . A 
even more dramatic reduction of the number of ONSs 
could follow from the recent investigation of ROSAT 
sources in dark clouds and high latitude molecular 
clouds by Danner (1996). For these fields predictions 
seem to be in excess by a factor ~ 10 — 100. 

Although the paucity of ONSs may be partly re- 



lated to the heating of the ISM by the source itself 
(Blaes, Warren, & Madau 1995), it could be also 
closely connected to the properties of NSs at birth 
and/or to the long term evolution of their physical 
parameters. Present estimates on the observability of 
ONSs rely on a number of physical hypotheses con- 
cerning the velocity distribution of pulsars at birth 
and the evolution of the magnetic field (e.g. Treves, 
& Colpi 1991; Blaes, & Madau 1993; Colpi, Cam- 
pana, & Treves 1993; Zane et al. 1995). The velocity 
distribution of ONSs was derived from the evolution, 
over the lifetime of the Galaxy, of the pulsar velocity 
distribution of Narayan, & Ostriker (1990, see also 
Paczyhski 1990) while a relic magnetic field <; 10 10 G 
was assumed a priori. 

Disappointingly, both the velocity distribution of 
NSs and the long-term evolution of the magnetic field 
are still affected by large uncertainties which make 
them highly controversial issues. The actual number 
of detectable ONSs depends crucially on the number 
of low velocity neutron stars, which are those accret- 
ing at the highest rate. This, in turn, is directly re- 
lated to the NS velocity distribution at birth. Lyne, 
& Lorimcr (1994) and Hansen, & Phinney (1998) de- 
rived a pulsar velocity distribution containing fewer 
slow objects than in Narayan, & Ostriker (1990), but 
this result is in contrast with other recent investiga- 
tions by Hartman (1997) and Hartman et al. (1997), 
which indicate that the number of pulsars with v <; 40 
kms -1 is the same as in Narayan & Ostriker or even 
higher. 

Observations show that neutron stars with mag- 
netic fields <C 10 12 G, possibly 10 9 G, are definitely 
present in LMXBs and in millisecond pulsars, suggest- 
ing that the magnetic field have decayed in these sys- 
tems. Instead, no firm conclusions have been reached 
yet for isolated objects (see e.g. Srinivasan 1997; 
Wang 1997). Theoretical results are far from being 
univocal and predict either exponcntial/power-law 
field decay (Ostriker & Gunn 1969; Sang & Chan- 
mugam 1987; Goldreich & Reisenegger 1992; Urpin, 
Chanmugam & Sang 1994; Miri 1996; Urpin, & Mus- 
limov 1992) or little or no decay at all within the age 
of the Galaxy (Romani 1990; Srinivasan et al. 1990; 
Goldreich & Reisenegger 1992; Pethick, & Sahrling 
1995; see also Lamb 1991 for a revue). Statistical 
analyses based on observations of isolated radio pul- 
sars give equally ambiguous results (Narayan & Os- 
triker 1990; Sang & Chanmugam 1990; Bhattacharya 
et al. 1992), owing in part to the difficulty of treat - 
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ing selection effects (Lamb 1992). At the present, 
there is no clear evidence of field decay during the 
pulsar phase, but this does not preclude the possi- 
bility of field decay over longer timescales. Different 
approaches to pulsar statistics led, independently, to 
the conclusion that, if the magnetic field decays, then 
it probably does so over a timescale ~ 100 Myr, well 
above the characteristic pulsar lifetime (Srinivasan 
1997, Hartman et al. 1997; Lorimer, Bailes, & Harri- 
son 1997). 

Very recently, however, Wang (1997) and Konenkov, 
& Popov (1997) noticed that a substantial decay of 
the magnetic field should have occured precisely in 
the ONS candidate RX J0720.4-3125. This object 
may have nevertheless a different origin resulting from 
common envelope evolution in a binary system. In 
the last hypothesis the decay of the -B-field may have 
been induced by accretion (see e.g. Shibazaki et al. 
1989; Romani 1990). 

Keeping in mind these uncertainties, the case of 
RX J0720. 4-3125 may provide, for the first time, evi- 
dence for a decay of the magnetic field in aging, iso- 
lated NSs. If evolution leaves a relic field ~ 10 8 — 10 9 
G a large fraction of the entire NS population is ex- 
pected to accrete the ISM, since major spin-down by 
propeller occured over the lifetime of the stars, as 
pointed out by Blaes, & Madau (1993). Actually, 
previous investigations assumed that nearly all ONSs 
accrete from the ISM (see again Treves, & Colpi 1991; 
Blaes, & Madau 1993). At the low rates typical of 
ONSs embedded in the ISM, however, the conditions 
under which accretion is possible are rather stringent, 
since the inflowing gas may or may not penetrate in- 
side the accretion and the Alfven radii according to 
the star velocity, spin period and magnetic field. As 
a consequence the number of ONSs accreting today 
may depend sensibly on the distribution of the stellar 
parameters at birth but, even more, on the long-term 
variation of the spin and magnetic field. 

The aim of this investigation is to explore the role 
of the magnetic field evolution in affecting theoret- 
ical predictions for ONSs observability. In order to 
do this, we relax the assumption of a constant, relic 
magnetic field and consider a model in which spin 
evolution causes the core field to migrate to the crust 
where dissipation processes (such as electron-phonon 
scatterings and scatterings on impurities) drive the 
field decay (see e.g. Srinivasan et al. 1990; Urpin, 
& Muslimov 1992). Several evolutionary tracks have 
been computed and the properties of the ONS candi- 



date RX J0720.4-3125 are easily recovered. We show 
that the paucity of ONSs detections can be accounted 
for if the magnetic field decays on a characteristic 
timescale ~ 10 s — 10 9 yr. The approach presented 
here suggests that, in close analogy with what is done 
for pulsars, the study of ONSs statistics may reveal 
a promising probe for constraining the spin and mag- 
netic field evolution of neutron stars. 

2. Spin and Magnetic Field Evolution 

NSs are widely believed to be borne with high mag- 
netic fields, B « 10 12 G, and very short rotational 
periods P « 10~ 2 s. During the early phase of the 
evolution, magnetic dipole losses spin down the star 
at a rate which, in vacuo, is 

^lO-^P^syr- 1 , (1) 

where B 12 is the surface dipole field in units of 10 12 
G. The NS mass and radius have been assumed to be 
M* = IMq and R* = 10 km throughout. Actually, 
isolated NSs are surrounded by the diffuse interstellar 
medium, so the pulsar relativistic wind continues un- 
til the ram pressure of the infalling material exceeds 
the outflowing momentum flux at the accretion radius 
(e.g. Blaes, & Madau 1993; Treves, Colpi, & Lipunov 
1993 and references therein). This occurs when the 
period exceeds a critical value 

Ppror-Wn-^vKfBlfs, (2) 

where v 10 = {V? + c 2 10 ) 1/2 , V w is the ONS velocity 
in units of 10 km s _1 , c s ~ 10 kms~ is the ISM 
sound speed and n the ISM number density in cm~ 3 . 
Here and in the following the Bondi formula for the 
accretion rate 

M ~ lO'V^gs" 1 (3) 

is used. For P > P prop the incoming gas can penetrate 
down to the Alfven radius 

R A ~ 1.7 x 10 10 i^ 2 /7 n-^vli 7 cm . (4) 

The interaction of matter with the magnetosphcre 
prevents, however, accretion to go any further, be- 
cause, for typical ONSs parameters, the Alfven radius 
exceeds the corotation radius R co = 1.5 x 10 8 P 2 / 3 cm. 
Since the flow is forced into supcrkeplerian rotation, 
angular momentum will be extracted and spin-down 
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proceeds at a rate (the propeller phase, Illarionov, & 
Sunyaev 1975) 



dP_ 
~dt 



10- 8 n 9 / 13 «- 27/13 Si2 /13 ^ 21/13 s yr' 1 . (5) 



Only when the period exceeds 



P a ~ 2.5 x 10 J B^ 2 ' n 



3 R 6/7 -3/4 l/7 t 



(6) 



steady accretion onto the star surface can set in. 

^From equations ((fl), @, (||) and (||) both the crit- 
ical periods, P pr . op and P a , as well as the spin-down 
rates, depend on the magnetic field. Hence the overall 
picture of the NS evolution would change significantly 
if the magnetic field is itself evolving in time. In this 
paper we consider a simple decay model, in which the 
magnetic field is anchored to both the superfluid core 
and the crust, so that its evolution is primarily con- 
trolled by the star rotation (see e.g. Ding, Cheng, 
& Chau 1993; Miri, & Bhattacharya 1994). A differ- 
ent model, in which the field is confined in the star 
crust will be briefly discussed in section 5. In the 
present picture, the neutron vortices, carrying the an- 
gular momentum, are pinned to the proton fluxoids. 
The radial drift of the vortices, induced by spin-down, 
causes the magnetic flux to migrate to the crust where 
ohmic dissipation occurs (Srinivasan 1990; Srinivasan 
1997). The physics of the interaction between vortices 
and fluxoids is still uncertain, but, as shown by Ding, 
Cheng, & Chau (1993) (see also Miri, & Bhattacharya 
1994), the decay rate of the core field B c turns out to 
be approximately equal to the spin-down rate 



dB c 
dt 



dP 
~dt 



(7) 



yielding B c = B c (0)P(0)/P. Ohmic dissipation takes 
place in the crust inducing the decay of currents on 
a characteristic timescale r c (see again Ding, Cheng, 
& Chau 1993; Miri, & Bhattacharya 1994). This re- 
sults in the evolution equation for the surface mag- 
netic field B 

dB _ B-B c 

The crustal timescale r c depends on the mechanism 
responsible for ohmic dissipation and is, in general, a 
function of time. However, in the spirit of the simple 
model used here and because of the still poorly known 
details of the microphysics, we found more convenient 
to view t c as a constant which we vary to bracket 



uncertainties. From equation QSJ) one can see that the 
core and surface fields evolve on the same timescale if 
r c is shorter than the spin-down time (controlling the 
evolution of B c ). In the opposite case, the core field 
decays faster while the surface field relaxes to B c on 
the scale r c . 

3. Do Aging NSs Accrete Steadily ? 

In this section we follow the evolution of the spin 
and the surface field integrating numerically equa- 
tions (|l|), (H) and (||). At the onset of evolution, 
at t = 10 3 yr, the neutron star is assumed to spin 
rapidly with a period P(0) = 0.01 s, while the ini- 
tial value of the surface field is taken to be consistent 
with the narrow gaussian distribution (log -Bo = 12.4, 
a = 0.32) inferred from pulsar observations (see 
e.g. Bhattacharya et al. 1992). Accordingly, we 
restrict our analysis to values of B(0) in the range 
8 x 10 11 < B(0) < 7 x 10 12 G. The initial core field 
B c (0) is taken to be always equal to B(0)/2 (Ding, 
Cheng, & Chau 1993). 

The evolution is divided into three phases (see e.g. 
Blaes, & Madau 1993; Treves, Colpi, & Lipunov 1993; 
Wang 1997). As long as P < P pr0 p, rotational en- 
ergy losses are due to magnetic dipole radiation and 
the behaviour of B and P is governed by equations 
and (|). For P prop < P < P a , equation (Q) is 
replaced with ([s]) to account for spin-down by pro- 
peller. Since the latter depends on the accretion rate, 
we vary the star velocity V and keep the ISM density 
constant (n — 1cm -3 ), in exploring different evolu- 
tionary tracks. The numerical integration is halted 
either when the period P becomes equal to P a (cal- 
culated with the current value of B) after a time t a 
or when t = to = 10 10 yr, the age of the Galaxy. As 
discussed by Wang (1997), during the accretion phase 
the mass loading of the field lines could become im- 
portant in spinning down the star, because it increases 
the moment of inertia of the star/magnetosphere sys- 
tem. However, in the case of RX J0720. 4-3125 and 
exponential decay, it only produces a small variation 
(from 7.86 s to 8.39 s) in the value of P. Since here 
we focus on the average properties of the ONS pop- 
ulation, we neglect this further source of spin-down, 
assuming that P (and, accordingly, B c ) is a constant 
for t > t a . The surface magnetic field, however, con- 
tinues to decay and equation (0) gives 



B (t) = (B c ) a + [B* - [B c ) a ] exp 



t-t* 



(9) 
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where (B c ) a = B c (0)P(0)/P a . In the following we 
consider different values of the decay time t c and dis- 
cuss the ensuing evolutionary scenarios. 

3.1. Very Fast and No Decay 

Although statistical studies of isolated radio pul- 
sars indicate that little or no field decay occurs over 
the pulsar lifetime (Bhattacharya et al. 1992; Lorimer, 
Bailes, & Harrison 1997), values of r c ~ 10 7 yr can 
not be completely ruled out. Anyhow, we find that 
rapid crustal field decay inhibits angular momentum 
losses to such an extent that virtually no ONS would 
be accreting today, all of them would either be dead 
pulsars or stay forever in the propeller phase. Quite 
the opposite happens if r c = 10 10 yr, a value large 
enough to mimic the case of no field decay in the pop- 
ulation of ONSs. In this case most neutron stars enter 
the accretion stage remaining highly magnetized and 
rotating very slowly, with periods up to few days. 

3.2. Decay on Intermediate Timescales 

We start our discussion examining in some detail 
the evolution of P and B for r c = 10 8 yr. Figure la, 
b illustrates the history of a low velocity neutron star 
with V w = 2 and log B(0) = 12 A G. The endpoint 
of the evolution has properties similar to those of the 
neutron star in RX J0720.4-3125 (Wang 1997). As it 
can be seen, during the lifetime of the Galaxy the star 
goes through all the three phases, pulsar, propeller 
and accretor, and may be visible at present as an 
accreting source. The situation is entirely different 
for a NS with higher speed, Vio = 30. In this case, 
after the initial loss of rotational energy by magnetic 
dipole radiation, the star enters the propeller phase, 
but has no time to spin down sufficiently to become 
an accretion-powered X-ray source because M is too 
small to make braking effective. This means that NSs 
borne with high kick velocities may undergo only the 
pulsar and propeller phases over their entire lifetime. 

Neutron stars which become accretors arrive at the 
end of the propeller phase with a period P a and a 
surface field B a which depend on the accretion rate, 
and hence on the star velocity, once n, B(0) and r c are 
fixed. The behaviour of B a and P a as functions of V 
is shown in figures 2a, b where each curve corresponds 
to a different value of B(Q). In figure 2a triangles and 
squares mark the loci where t a = 10 9 and 10 10 yr, 
respectively. 

As it can be seen, for each values of B(0) there 



exists a limiting velocity, V cr , above which the time 
required to reach P a exceeds the age of the Galaxy. 
Clearly V cr depends on B(0): NSs with large fields 
at birth can reach the accretion phase even if they 
move with comparatively large speed, while for weak 
initial fields rapidly moving objects remain in the pro- 
peller stage. We notice (see figure 2a) that as V in- 
creases toward V cr substantial field decay has occured 
before the accretion phase begins, and B a approaches 
~ 10 9 G, irrespective of J5(0). For V ~ V cr , B a repre- 
sents the relic field, since t c <C t a ~ to- The distribu- 
tion of the periods is rather wide and ranges from few 
seconds to several minutes (see figure 2b) . For low val- 
ues of B(0) (bottom curve of figure 2b) P a monoton- 
ically increases with V, while at larger B(0), P a (V) 
exhibits a maximum. This is a not in contrast with 
equation (||) because P a depends also on B a which, 
in turn, is a function of V. 

In our model the periods at t = t a coincide with 
those at t = tc, while this is not the case for the 
surface field because of further ohmic dissipation for 
t a < t < t(j. The distribution of the surface field at 
present, B now = B(to), is plotted together with the 
core field in figure 3a. As it appears from figure 3a, 
decay is fast enough to freeze the surface fields to the 
core field at t = t a . ONSs with spin periods ^5 — 10 
minutes and low magnetic field (B ~ 10 s G) can be 
easily accommodated by the model for r c = 10 s yr 
(see figure 3b). 

Ohmic dissipation of the crustal field on a scale 
~ 10 9 yr speeds up the evolution to the final state 
of accretor. Both losses by magnetic dipole braking 
and propeller slow down the NS effectively. Conse- 
quently, V cr takes higher values and also NSs with 
large enough velocities can enter the accretion stage. 
Spin periods vary between few seconds and several 
hours and, despite the longer crustal decay time, 
B now may be lower than for r c = 10 8 yr. This is 
because the larger spin-down rate produces a signifi- 
cant decay of the core field and there is enough time 
after the onset of accretion for the surface field to 
approach the core field. Typical values of B now are 
between 10 8 and 10 9 G (see figure 4a). We notice 
that for r c = 10 9 yr the period of RX J0720.4-3125 
can not be reproduced along any of the evolutionary 
tracks (see figure 4b). 
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4. More On ONSs Detectability 

Within our spin-induced field decay model, we now 
estimate the fraction of NSs (relative to the total 
number) that are presently in the propeller, N prop , 
and in the accretor phase, N a . To this end, we recall 
that only those stars having V > V cr (V < V cr ) are in 
the propeller (accretor) phase. N prop and N a depend 
on the velocity distribution of ONSs G(V,i), which 
changes during the evolution of the Galaxy. Starting 
from the pulsar distribution of Narayan, & Ostrikcr 
(1990), Blaes, & Madau (1993) and Zane et al. (1995) 
derived the present distribution function following the 
evolution of a large number of stars moving in the 
Galactic potential. The current cumulative velocity 
distribution may be conveniently approximated as 



G(V,t G ) 



1 + X r ' 



(10) 



where x — Vio/6.9, m = 3.3 and G(V,ta) gives the 
fraction of stars with V < V. The stars enter the 
propeller phase at t > 10 8 yr, so there have been 
time enough for their velocity distribution to evolve. 
In sections 2 and 3, we neglected this effect and as- 
sumed that V is constant in time. We judged this 
to be a minor flaw with respect to other sources of 
uncertainty present in our model, including the decay 
model itself. 

Table 1 contains N a and N prop for different r c and 
for log-B(O) = 12.4. Previous models relied on the 
hypothesis that all ONSs are currently accreting, i.e. 
N a = 1. Here we find that only for r c ~ 10 10 yr, 
most (~ 85%) of the NSs are now accreting. This fig- 
ure drops significantly to ~ 5% if we allow for crustal 
field decay on a timescale r c ~ 10 s yr. However, the 
fraction of potentially observable NSs is not N a , but 
depends on the characteristics of the detector. Here 
we refer to ROSAT, and the fourth column of table 
1 gives the fraction Nt of accreting NSs with emis- 
sion falling in the T-band (0.1-2.4 keV). The energy 
at which a blackbody spectrum peaks was taken as 
representative of the mean photon energy and the ef- 
fective temperature of radiation was estimated con- 
sidering polar cap emission (see figure 5). It is of 
interest to notice that, although the fraction of ac- 
creting NSs increases with r c , the relative number of 
sources in band is smaller for t c — 10 9 yr. The rea- 
son is that the surface field has decayed more (see 
figure 4a) than in the case r c ~ 10 8 yr, so the emit- 
ted spectra are softer. We find that, for t c ~ 10 8-9 
yr, Nt clusters around ~ 4%. Yet, this number can 



not be compared with previous figures on the num- 
ber of detectable ONSs (see Treves, & Colpi 1991; 
Zane et al. 1995) since they evaluated the number of 
sources with a flux above the sensitivity threshold of 
ROSAT. In order to obtain a meaningful comparison, 
we have computed the fraction N T d of ONSs in the 
T-band, assuming a magnetic field of 10 9 G, n = 1 
cm -3 and blackbody polar cap emission. N T d — 0.21 
gives precisely the number of sources in the T-band 
in the spirit of previous models and R = Nt/N^ 
is the ratio between the present and past predicted 
numbers. The values of R are given in Table 1 and 
indicate that the expected number is reduced by a 
factor ~ 4 for r c = 10 8 yr. R is even smaller for 
t c = 10 9 yr owing to the more severe decay of the 
core field. For t c ~ 10 7 yr, R is as low as 0.03; by 
contrast, R exceeds unity for r c = 10 10 yr. We note 
that N prop and N a are only very weakly dependent 
on the initial value of the period: models with larger 
values of P(0) give essentially the same figures. 

5. Discussion and Conclusions 

In this paper we have shown that the expected 
number of ONSs observable with ROSAT can be sen- 
sibly reduced if the star magnetic field decays on 
timescales ~ 10 8 — 10 9 yr. ^ This is primarily due to 
the lower fraction of NSs which are in the accretion 
phase now, but also to a softening of the spectra which 
may bring the peak of the emission outside ROSAT 
T-band if r c = 10 9 yr. It follows that, if the paucity 
of detected ONSs is to be explained in terms of a 
global reduction in the number of accretion-powered 
sources, decay times r c <J 10 8 yr, should be invoked. 
For intermediate values of r c , ~ 10 9 yr, the fraction of 
accreting ONSs is ~ 50 — 60%, too large to match ob- 
servational limits, but only ~ 15% of the sources falls 
in the ROSAT T-band, yielding again an acceptable 
number of detectable objects. This last figure should 
be regarded as indicative, being based on the assump- 
tion that all sources emit a blackbody spectrum from 
the polar caps. If the decay time is ~ 10 10 yr or larger, 
spin-down by propeller is so effective that nearly all 
ONSs are accreting today. In this case the surface 
magnetic field is close to the initial value, B ~ 10 12 G, 
and deviations from a pure blackbody spectrum are 
to be expected. Unfortunately, no detailed radiative 
transfer calculations are available for accretion atmo- 



2 Soon after this paper was submitted, we became aware that a 
similar conclusion was reached by Livio, Xu, & Frank (1998). 
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spheres onto magnetized NSs, at least to our present 
knowledge. Some hints may be, nevertheless, derived 
from the analysis of radiative transfer in cooling neu- 
tron stars atmospheres (Shibanov et al. 1992), since 
the input physics of the two model is quite similar. 
The emerging picture indicates that magnetic effects 
tend, somehow, to balance the hardening found by 
Zampieri et al. (1995) in the case of unmagnetized 
accretion atmospheres, so that, for B ~ 10 12 G, spec- 
tra are, globally, not too different from a blackbody 
emitted by the polar caps. With this proviso, figures 
in table 1 tell us that all accreting sources should be 
indeed within ROSAT T-bandpass and show that de- 
cay times t c ;> 10 10 yr provide a fraction of potentially 
detectable ONSs equal to that derived in previous in- 
vestigations. 

Our main conclusion is that models of spin-induced 
magnetic field decay with a crustal dissipation time 
shorter than the age of the Galaxy are able to bring 
theoretical predictions of ONSs detectability closer to 
present observational limits. It is interesting to note 
that ohmic diffusion times between 10 8 5 and 10 9 yr 
can also explain the low values of B observed in mil- 
lisecond pulsars, which experienced spin induced field 
decay during binary evolution in wide low-mass sys- 
tems (Miri, & Bhattacharya 1994). 

Although our present results seem indeed to point 
to a decay of magnetic field in ONSs, a firmer con- 
clusion may be reached once the number of sources 
actually present in ROSAT fields is computed evolv- 
ing the ONSs distribution function together with the 
spin and the magnetic field. Two other issues deserve 
further discussion. In computing the spin-down rate 
in the propeller phase the original expression of II- 
larionov, & Sunyaev (1975) was used. It should be 
noticed that the physics governing the propeller is 
rather complex and still far from a complete under- 
standing. It has been shown (see e.g. Ghosh 1995) 
that the characteristic time for angular momentum 
losses strongly depends on the details of the accretion 
process and can be longer than what is assumed here. 
If the propeller is less efficient in braking the star, the 
number of accreting ONSs reduces further. On the 
other hand numerical simulations (Wang, & Robert- 
son 1985) seem to indicate that propeller spin-down 
can go even faster that in equation (|^). 

The second point concerns the field decay model. 
If the field is confined only in the stellar crust (see e.g. 
Urpin, & Muslimov 1992) the emerging picture may 
change since the magnetic field is not coupled to the 



period. In the model of Urpin, & Muslimov t c is com- 
puted from the microphysics and varies in time, de- 
pending on the thermal evolution of the stellar crust. 
We have repeated the calculation using equation (||) 
with B c = and r c estimated from figure 5 of Urpin, 
& Muslimov, and found that the number of accreting 
ONSs again reduces. Even if the initial surface field 
is as large as 10 13 G only <; 1% of the NS population 
is the accretion phase now. On the opposite, an ex- 
ponential decay of the surface field with t c ^ 10 9 yr, 
as inferred by Pethick, & Sahrling (1995) in the case 
of a non-vanishing core field, do not seem to solve the 
problem of the excess of potentially visible ONSs. 

The conclusion that many ONSs should be in the 
propeller stage seem rather robust if the magnetic 
field decays over a characteristic r c ~ 10 9 yr, irre- 
spective of the details of the model. Is there a way to 
observe directly these objects ? The power associated 
to the propeller derives ultimately from the rotational 
energy, and therefore one can deduce that it should be 
smaller than that released by accretion of matter onto 
the star surface typically by a factor ~ 100 (see e.g. 
Treves, Colpi, & Lipunov 1993). The accreted mass, 
however, could accumulate at the Alfven radius un- 
til the pressure required to overcome the centrifugal 
barrier is built up, giving raise to recurrent episodes 
of accretion. If this is the case, ONSs in the pro- 
peller phase could be observed as transient sources 
(Treves, Colpi, & Lipunov 1993). More definite pre- 
dictions about the observational appearance of these 
sources, including spectral properties, luminosity, re- 
currence times, demand, however, for a more realistic 
modelling of the physics governing the propeller. 

We are grateful to Sergio Campana for construc- 
tive comments during the early stages of this work 
and to Piero Madau and Andrea Possenti for valu- 
able discussions. 
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Fig. 1. — a) The evolution of the surface field for 
Via = 2, log 5(0) = 12.4 G and r c = 10 s yr. b) 
The evolution of the period for the same choice of the 
parameters; the periods for the onset of the propeller 
(eq. (2| . dashed line) and of accretion (eq. || , dash- 
dotted line) are also shown. 

Fig. 2. — a) The surface field at the onset of 
accretion as a function of the star velocity for 
t c = 10 s yr, from bottom to top: B(0) = 
0.87, 1.37, 2.06, 3.06, 4.60, 7.28 x 10 12 G; triangles and 
squares mark the loci where t a = 10 9 and t a — 10 10 
yr respectively, b) Same for the spin period. 

Fig. 3. — a) The the surface (full lines) and the core 
field (dashed lines) at present vs. V (B c and B ac- 
tually coincide so the two curves are superimposed); 
details as in figure 2. b) The value of the surface field 
at present vs. the period; the vertical line marks the 
period of RX J0720.4-3125. 

Fig. 4. — Same as in figure 3 for r c = 10 9 yr. 

Fig. 5. — The effective temperature against the star 
velocity for log_B(0) = 12.4 and different values of 
r c ; the two horizontal lines mark the ROSAT T- 
bandpass. 
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Table 1 

Fraction of objects in the different phases* 



t c (yr) 


N a 

1 v prop 


N b 




R d 


10 8 


0.95 


0.05 


0.05 


0.23 


10 9 


0.37 


0.63 


0.03 


0.16 


10 10 


0.15 


0.85 


0.85 


4.10 



*Here log 5(0) = 12.4 
a Fraction of NSs in the propeller phase 
b Fraction of NSs in the accretion phase 
c Sourccs in the ROSAT T band 

d Ratio of sources in the T-band ac- 
cording to present and past models 
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